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The classic problem of the growth of density perturbations in an expanding Newtonian universe 
is revisited following the work of Bisnovatyi-Kogan and Zel'dovich. We propose a more general 
analytical approach: a system of free particles satisfying semi-degenerate Fermi-Dirac statistics 
on the background of an exact expanding solution is examined in the linear approximation. This 
differs from the corresponding work of Bisnovatyi-Kogan and Zel'dovich where classical particles 
fulfilling Maxwell-Boltzmann statistics were considered. The solutions of the Boltzmann equation 
are obtained by the method of characteristics. An expression for the damping time of a decaying 
solution is discussed and a zone in which free streaming is hampered is found, corresponding to 
wavelengths less than the Jeans one. In the evolution of the system, due to the decrease of the 
Jeans length, those perturbations may lead to gravitational collapse. At variance with current 
opinions, we deduce that perturbations with A > Aj Max/1-48 are able to generate structures and 
the lower limit for substructures mass is M = Mj max /(1.48) 3 « Mj max/3, where Mj max is 
the maximum value of the Jeans mass. 
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INTRODUCTION 

Bisnovatyi-Kogan and Zel'dovich (1970) |l| approached analytically the study of gravitational 
clustering of classical particles within the framework of a Newtonian universe. The general rel- 
ativistic approach is necessary only for perturbations with wavelengths A > Xnorizon (see e.g., 
Peebles (1993) ^)). They considered the evolution of density perturbations for particles obeying 
Maxwell-Boltzmann statistics. They found that for perturbation lengths A 3> \j the regime is the 
hydrodynamical one, while for A <C A./, they found a damping of density fluctuations instead of the 
traditional oscillatory behaviour for perturbations in a fluid (see e.g., Binney and Tremaine (1987) 
Q): they pointed out that short wavelengths perturbations are damped due to some purely kinetic 
mechanism. They studied only small and very large wavelengths with respect to the Jeans length 
and treated the Jeans length of the perturbation as the limiting length scale between the stable per- 
turbations and the unstable ones. The cosmological large-scale structure formation scenario in the 
work of Bisnovatyi-Kogan and Zeldovich (devoted to the study of a typical astrophysical collisionless 
system, like the edge of stellar clusters or clusters of galaxies) was the framework in which the self- 
gravity of collisionless particles (stars or galaxies) is in competition with the universal expansion to 
create a condensation of matter. This approach was not taken in order to study a top-down scenario 
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in which big perturbations have successive fragmentations, even if the Zel'dovich top-down scenario 
of pancakes was proposed during those years |5| . 

On the other hand Gilbert (1966) H, numerically solving an integral equation, has claimed 
that collisionless damping is similar to the Landau damping in plasma physics. Recall that in 
plasma physics the Landau damping is a phenomenon in which the energy of a definite wavelength 
perturbation is efficiently dissipated by collective motions of the system M. 

Bisnovatyi-Kogan and Novikov (1980) || and Bond, Efstathiou and Silk (1980) ^ have discussed 
the physics of neutrino masses in the expanding universe. When massive neutrinos are relativistic 
their Jeans mass Mj is larger than the horizon mass Mh, while when they become nonrelativistic 
the corresponding Jeans mass becomes smaller than the horizon mass, making possible collapse of 
masses Mh > M > Mj due to gravitational instability. They obtain the analytic formula for the red 
shift z n . r .( m v) x m v 2 at which the Jeans mass of the neutrinos enters the horizon, at this moment 
the whole Jeans mass is causally connected and its gravitational collapse starts (see fig. 2). The 
relation between mass and the dimensions of large scale structure can be recovered: scaling to the 
present radius of the universe Xj n.r.{rn v ) ~ 120Mpc( ™£ v ) ~ 1 . The neutrino is the first candidate 
for the dark matter, and leads naturally to a top-down scenario, owing to its low mass value. 

The work of Bond and Szalay (1983) §, based on a numerical approach, assessed the damping of 
the perturbations smaller than the Jeans length. They have followed a general relativistic approach 
based on the Gilbert (1966) method j|. This method contains the first order approximation of the 
equations of motion as in the Bisnovatyi-Kogan and Zel'dovich one. 

White, Davis and Frenk (1984) [H have shown that free streaming should only permit the onset 
of very large and hot clusters of gas (using the fraction of baryonic and therefore visible matter as 
a one-to-one tracer for the dark matter distribution) at KT w 20KeV and higher, too luminous in 
X-rays with respect to the observed ones at this epoch (the hottest known is the Perseus cluster at 
KT = 7KeV). 

The alternative cosmological model to a neutrino-dominated universe appeared in 1984 with the 
introduction of Cold Dark Matter by Primack and Blumenthal JTTJ , and the Hot Dark Matter models 
(including neutrinos) began to be somewhat neglected |l|. 

We have used the method of Bisnovatyi-Kogan and Zel'dovich in a detailed analytic study of 
the fragmentation for gravitational instability due to the Jeans mass of fermions, extending their 
approach with care to wavelengths on the order of the Jeans length. We find a more general 
expression for the damping time formula. 

We have evaluated the role of free streaming by studying the solution of the gravitational Poisson 
equation for a fermion density perturbation in the expanding universe. We have obtained a new lower 
limit for the free streaming length that is 1.48 times smaller than the Jeans length. The fundamental 
role for the process of structure formation appears to be linked to this new length, determined by 
the dynamics of the self-gravitating fermion system, and not to the traditional free-streaming length 
determined only by the kinematics of the particles. 

In section I, we present the mathematical tools for treating the growth of fermion density per- 
turbations in the expanding universe. In section II, we review the free-streaming effect, while in 
section III, we define the damping criteria applied to the decaying solution of section I. In section IV, 
the properties of the nondecaying solutions belonging to the instability shell are discussed. Finally 
in section V, we outline the conclusions. The appendix applies the method of characteristics to 
integrate the Boltzmann equation. 



I. SOLUTION OF THE LINEARIZED BOLTZMANN EQUATION 

A kinetic treatment is necessary to analyze the stability of a system of weakly interacting massive 
particles after their decoupling when the time between collisions is much greater than the character- 
istic hydrodynamic time (i.e., the crossing time t c = — , where I is the linear dimension of the system 
and v B is the velocity of sound or, better, the average sound velocity). We consider a self-gravitating 
system of particles in the background of an exact expanding solution and we examine its stability 
in the linear approximation. 

We describe the collisionless matter with a distribution function f(x, v, t) satisfying the kinetic 
Boltzmann equation 

and the Poisson equation 
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V 2 $ = AirGp, 

where G is the gravitational constant and the total density at x is 

P{%) = / f(x, v, t) dv. 



We study the development of perturbations in a Newtonian universe with critical density con- 
taining collisionless particles. 

It is known that for such an universe 

po = ^ ; $o = r°p^ 2 

(see e.g. Peebles (1993) Q). Solving the Boltzmann equation by the method of characteristics 
(Appendix) , we find that the comoving velocity is an integral of the motion, so that the solution fo 

of the unperturbed Boltzmann equation for the distribution function is an arbitrary function of the 

T • 2 2 - 1 

comoving velocity m = Vits — ^xtt 3 : 

fo = af(u). 

The constant a is determined by imposing the condition of self-consistency H: 

fod 3 v ' 



6-TvGt 2 



The distribution function describing the collisionless matter is f(x, v, t) such that / = /o + fi 
and the gravitational potential is $ = $o + $i (/i and "l?i are small perturbations) satisfying the 
kinetic Boltzmann equation. Following Bisnovatyi-Kogan and Zel'dovich jij we seek the solutions 
in the form 

$i = e* s >(t) , ft = e*«/(t), 

where ^ = x/tz is the comoving coordinate and k = (fc,0,0). The solution of the linearized kinetic 
Boltzmann equation is 

/r= T^§l<A(0) = 0. 
j ox' ov' 

Substituing the expression (Q) into (^), using the relations (A13), (A14) for x\ and v[ in the 
Appendix, and integrating by parts, we have 

Here ui is the component of the comoving velocity u along the direction of the comoving wave 
number k. The Poisson equation, using the expression (Q) becomes 

t3 J 

Then integrating Equation (Eh over the velocities and substituting the result into eq.(|lo[), we obtain 
the following integral equation for <f>: 

16tt 2 G rt 



t*(j>(t)-\ — / (/)(t )dt du ufa sin(3fcitr) = 0, 

k Jo Jo 

where u = Ui is the comoving velocity, r = t~s — t' - 3, and k is the comoving wave number. We 
are considering the universe with critical density filled by fermions, decoupling in the relativistic 
regime, so that the unperturbed semi-degenerate Fermi-Dirac distribution function fo is 
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8niG<u>3g(O e ^>-i + 1 
where 



^ 3k B T v (t) 

< u >= , 

mc 

and £ = -jjAjt is the degeneracy parameter, ks is the Boltzmann constant and 11 is the chemical 

potential. The normalization function is g(£) = |£ 3 +2£r(2)£ + 2 — e_s ^> where Cfl( n ) — 

^ s t ne Riemann function of index n. In equation ( |l3|) the relation between energy (here 
at equipartition) and momentum for those fermions is in the relativistic regime E — cp because we 
assume that they are decoupled from the radiation and the matter while relativistic. The fact that 
the particles are relativistic at decoupling only depends on the weak interactions and not on the 
fermion mass jl7|. For the Liouville theorem extended to the expanding universe fTf^, this relation 
holds at all times. In formula (^|) the comoving velocity u = iti along one direction is multiplied 
by three under the hypothesis of spatial isotropy. 

Inserting the distribution function ( |l^ ) in the integral equation (pi]), we obtain 
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where a = k<u>,y = and < u >=< v > 2 t 3 . Defining the Jeans length as 

2 1 , 7T . i 

A. 7 =<t; 2 >= — )« 
Gp 

at a given time t = to, using also p = g^gp , we can write 

2 it 1 i 
a = k < u >= — (— )t 3 = FatS 



where a = jj, and F = ^ « 1.18. Defining 



and 
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and 7] = x' — x, the dimensionless integral equation is 

ysin(a77y) 
e«"« + 1 



with 



r(a,£,t7) = / dy 



i) 



A numerical analysis shows that for fixed a and £, the function T increases and decreases rapidly 
with rj, having a maximum at rjo. 

The integral in the equation (|l9|) can be calculated by the method of steepest descents |^3j . We 
approximate T(a, £, rj) w e^iiO so that 

T? ( a; + 7 7) 4e (a; + r ? o) 4 |aft"('?o,0l i ' 

It is important to note that the function T(a,^,rj) is not a true function of three variables, but 
only a function of two variables: r(ar), £). Thus with the variation of the parameter a, the T 
function undergoes a conformal transformation along the rj axis: a contraction as long as a grows, 
and a dilation if it decays. The maximum of the T function at rjo for a given value of a is equivalent 
in the two- variable framework to a unique maximum at rj* = arjo = const (see fig. 1). So in the 
approximation of the T function by an exponential function in the steepest descent method, there 
are not distinct constraints for small or large values of Q, while Bisnovatyi-Kogan and Zel'dovich 
examine these two regimes separately. 

The integral equation ( Jig| ) becomes 



with 



ah"(r] ) 



f °° J„. V 3 sin(arjgy) 



Jo a V e »-4+i 



since ah(r/, £) = In r(a, 7;, £). In this approximation, as in the Bisnovatyi-Kogan and Zel'dovich work, 
we can Taylor expand the function (j>(x + 7/0) and (x + ryo) -4 at 7/0 = 0, retaining the terms of the 
first order in r/o, leading finally to the differential equation 

r) 4> (x) = 0, 



where 



with the solution 



B = V2^ 



a 2 g(£) x 



or in terms of time t, using 



and 



I8B770 



</»(<) = Cexp(l), 



7 =I8B^ F t0 - 



II. THE FREE STREAMING EFFECT 



We want to study whether the density perturbations smaller than the Jeans length will lead to 
substructures lighter than the Jeans mass, or will be dissipated by the collisionless damping. 

The free streaming kinematic scale A/ s (see fig. 2) is the scale that a particle with some mean 
velocity has traveled up to time t: 

X fs = / dt'a(t')v(t') 
Jo 

(where the scale factor of the universe a(t') oc (z + 1 

The horizon size Xnorizon is the distance traveled by a photon emitted at t = 0(see fig. 2). 

The Jeans scale Aj contains a mass of particles whose self gravity is greater than the kinetic 
pressure (see fig. 2). 

In the extremely relativistic regime for particles (see the left side of fig. 2) on scales larger than 
\f s ~ Xnorizon, the perturbations can neither decay nor grow p2| , while on scales A2 smaller than 
\f s and smaller than the Jeans mass scale Aj, the perturbation decays quickly due to directional 
dispersion. As soon as the particles become nonrelativistic (see the right side of fig. 2) , the pertur- 
bations should decay slowly due to velocity dispersion. 

When the Jeans length becomes smaller than the horizon length, the perturbations greater in 
scale than Aj and A/ 3 increase (see Ai on the right side of fig. 2 when it crosses the Jeans length 
curve) . 

Perturbations started after z n . r . with wavelengths A3 that are greater than A/ s (calculated inte- 
grating from to > t„. r .) will survive and collapse when crossing the Aj curve, while the ones like 
A4 will be dissipated by the free streaming. 

The current opinion Jio| , ^3| is that the free streaming scale is the lower limit for surviving pertur- 
bations that may lead to subsequent gravitational collapse. Below this scale, it is generally assumed 
that all perturbations are dissipated. In what follows we present an alternative approach to this 
problem. 



III. DAMPING CRITERIA 



The characteristic time for damping starting from the time to is the time interval T in which the 
density perturbation is reduced to | of its value at the time to. From the expression (E7J) one finds 



with 7' = and 

' In 



1' 



, = g(QF 3 k 2 

' 18B{Orio ( k.j' ' 



If the parameter 7' is near 1, the damping time can be very long and the density perturbations of 
the corresponding wave number k survive the collisionless damping process. 

The expression (^) is more general than the characteristic time of damping considered by 
Bisnovatyi-Kogan and Zel'dovich (1970) (equation 20). 

Recalling that ar]o = rf (Vo makes the numerator of B in ( pB| ) a maximum for a given value of a) 
and that a = where k j is the comoving wave number corresponding to the Jeans length Jl4| , we 
can derive the lowest value of the wave number k for which the equation (^) has negative values, 
corresponding to undamped solutions 



* _ 3 W^ L25 = 



k.r V 9ti)F3 

In the expression Q F w 1.18, w 1.80 for £ = 0, B w 0.60 and rf « 0.54, we have that 
fluctuations with ai < 1.25 are not damped. 

Another criterion defining the undamped perturbations is evident from fig. 3, in which we have 
represented the temporal evolution of the relative amplitude of the fermionic density perturbations 
j^, where 8 oc a(t)(f>(t) (see e.g. Peebles (1993) ||) and S = is the density contrast at time to- 
Setting 5 = leads to 

f = (^)exp[0.51a 3 ((f)-i-l)]. 
00 ao ao 

For each curve we have calculated the abscissa = (^j-) 3 , the expansion parameter of the 
universe) of the minimum point 

(-) « 0.83 a 2 . 

\ do / YfiAm 

For a curve attaining the value 
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we assume that the perturbation is undergoing a significant damping; numerically solving the 
trascendental equation ([35]) we obtain the value 

w 1.71 = q 2 ; 

from this second criterion the fluctuations with 02 < 1.71 are not damped. 
Summarizing the results found using these two criteria: 

A perturbation in the potential or in the density contrast function decays exponentially with time 
for short wavelengths. 

Furthermore in a given range of wavelengths near the classical Jeans length, this decaying slows 
radically and the perturbations stabilize. 

Due to the weak self-gravity the damping of the perturbations is faster for small linear dimensions 
(wavelengths). This behaviour is shown in the fig. 3 for high values of a. 

The physical mechanism responsible for the damping of the perturbations, in our model, is the 
free streaming effect. It is to be intended as a kinematical effect, whose amount depends on the 
statistics of the particles. 

Self-gravity is the competing mechanism for the growth of perturbations, which is a collective 
effect. 

The dynamical free-streaming scale results from these competing mechanisms and it also defines 
a length that is the lower limit for density perturbations capable of forming structures. 
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IV. THE PROPERTIES OF THE INSTABILITY SHELL AND PHYSICAL REMARKS ABOUT THE 

SOLUTIONS 



In figure 4 we describe the processes of growth, stability and decay of fermion density perturba- 
tions. The figure is divided into two parts using a vertical line that corresponds to the transition 
from the ultra-relativistic to the non-relativistic regime. 

On the left part of the figure we consider perturbations in the ultra relativistic regime. We select 
a perturbation of comoving length A2 which evolves along an horizontal line: for A2 > Xhorizon 
the perturbation is stable and has no damping since the different parts of the perturbation are not 
causally connected. For A2 < ^horizon < Xj eans the self gravity starts to act and the perturbation 
is quickly damped by free streaming. 

On the right part of figure 4 we consider the perturbations in the nonrelativistic regime. 

From the two damping criteria we have found that some perturbations of wavelengths smaller 
than both the Jeans length and the free-streaming length are stable and undamped. Selecting a 
perturbation with an intermediate value of = a between 1.25 (found with the first criterion) and 
1.71 (found with the second), we see that significant damping occurs only for ^- = a > 1.48. 

This is the definition of the instability shell, in which the perturbations with a < 1.48 are not 
damped and may lead to subsequent collapse due togravitational instability. 

So we can say that the free-streaming equation ( p9| ) describes the " kinematical" free-streaming, 
for which self gravity has no effect. The crossing of the perturbation length Ai (and A3) with the 
kinematical free-streaming length curves (that are extended as gray lines into the instability shell 
for showing the effect), in the right-hand region, does not imply the damping of the perturbation, 
because this length belongs to the instability shell. 

Our solutions satisfy the Poisson equation requiring a causal connection between particles, so 
they are not fully extendible to the relativistic regime, for which the Jeans scale is larger than the 
horizon. 

Meanwhile it is reasonable to assume a continuity of our solutions for the instability shell towards 
the relativistic regime. The instability shell becomes the region between the horizon and kinematical 
free-streaming lengths, and its lower limit goes continously from the kinematical free-streaming 
length to x \ ffis . 

Consequently we can deduce that perturbations with Ai > Xj Maa;/l-48 are able to generate 
structures. The lower limit for substructures (with respect to the maximum Jeans length) mass is 
M = Mj max /(lA8) 3 . 

Moreover for a starting perturbation of length A3 while the particles are nonrelativistic, the 
kinematical free streaming is not efficient, A3 is greater than the lower limit of the instability shell, 
and this substructure being stable, it can start the gravitational collapse when Xj becomes equal to 
A3. For A4 < A3 the perturbation is dissipated as in the classical case (see fig. 2). In the bottom 
right-hand part of figure 4 (for A3 and A4) one again finds on smaller scales and in the nonrelativistic 
epoch the same behaviour of the dynamical free-streaming length regulated by the instability shell 
for the primordial fluctuations Ai and A2. 

V. CONCLUSIONS 

The discovery of the instability shell, which depends on the statistical properties of the particles 
involved in the clustering phenomena, has clarified the role of free-streaming and furthermore leads 
to a better definition of it. We have introduced the concept of kinematical free streaming for the 
traditional treatment, which neglects self gravity and microscopical quantum statistics, while we 
introduce the new concept of dynamical free-streaming which takes into account these effects. 

Collisionless damping is a consequence of a purely kinematical effect (due to directional and 
velocity dispersion) of free streaming |l5) which does not act like a collective effect which the 
comparison with Landau damping (even the gravitational one Jl^| ) suggests. Free streaming is 
efficient only in the first stages of the damping process; after those initial transient stages the 
dynamics of the whole perturbation acts in the opposite direction, like a true collective effect. 
The successive damping inefficiency is a consequence of the microscopic statistics adopted for the 
particles: at the same temperature the fermions' velocity dispersion is less spread towards high 
values with respect to classical Maxwellian particles, thus the free streaming is less efficient for 
fermions in comparison with classical particles, as we have found. 
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The dynamical free-streaming length should refer to the instability shell's properties. Indeed, for 
semidegenerate Fermi-Dirac particles, primordial perturbations with wavelengths A > Xj 1 must 
survive damping process, this inequality determines the upper limit for the dynamical free-streaming 
length. 

We have identified an intermediate phase, with respect to the previous classical frameworks, 
in which the perturbation (within the horizon) neither grows nor decays until the corresponding 
Jeans mass is reached; at that point the collapse starts. This property allows structures smaller 
than the maximum Jeans mass to survive and makes possible their fragmentation into smaller 
substructures that can detach from the Hubble flow. 

Furthermore for an induced clump (for instance, hot gas from active galactic nuclei explosions 
or generated during processes of galaxy and cluster formation ]2^]) of lengths A > Aj/1.48, while 
the particles are nonrelativistic, the perturbation belongs to the instability shell. It will lead a 
substructure over a relatively rapid time scale, via the Jeans instability. 

Moreover the critical ratio a = 1.48 so obtained is independent of the scale of the perturbations 
considered, so this property is scale invariant and can be valid even after more fragmentations and 
may even lead to local fractal- like distribution of matter pi) . 



APPENDIX: 



The method of characteristics is applied to solve a quasi-linear and homogeneous equation of the 
form 

_, . du , .du . du 

in which the functions P(x, y, z), Q(x, y, z),R(x, y, z) are continuous and not simultaneously vanish- 
ing. Let us consider the continous vector field 



F = P{x, y, z)i + Q(x, y, z)j + R(x, y, z)k 

with i,j,k unit vectors along the coordinate axes. The characteristics of ( |Alj ) are integral curves of 
this vector field, defined by integrating the following equations 



dx 



dy 



dz 



P(x,y,z) Q(x,y,z) R(x,y,z)' 

Two independent integrals of these equations ^i(x,y,z) = Ci and $2(2;, y, z) = C2 represent 
surfaces whose intersection curves are these characteristics. The general solution of (Al) can then 
be expressed in the form u = $(Ci, C2). 

The characteristics of the Boltzmann equation jl[ 

21 + ^21 ^21^1 = 

dt dx dx dv 



are determined by the system: 



dt 



dxi 



dvi 



from which we obtain 



and 



xt = J Vidt, 



dvj 
dt 



2xj_ 
"9 t 2 ' 



Substituting eq.(<Y5) in eq.(A7) leads to 



and differentiating this with respect to the time t we obtain the second order differential equation 

dvi 2 



,2d 2 Vi .. , 

1 -dF + 2t Ht + 9 v ' =0 - 



If t — e z we have 



The solution is 



d Vi : dv % 2 



u« = Cut 3 - C 2 it~ 



Substituting eq. (A7) in equation (A5) and integrating, we have 



Combining the two expression (A.7) and (A8), we obtain the six integrals of the system (A4) 

Cii = SVifs - Xi t~i, 



C 2i =Vit 3 - -Xit 3. 



An arbitrary function f{C\i,Cii) of these integra ls is a solution of the transport equation (A3). 
The integral of the motion Cn given by eq. (A10) coincides with the comoving velocity u 



Solving the integrals (|A9|) and (A10) with respect to x' and v' we have 



from which we obtain 



and 



2 2 _i ,2 2 , _i 

Vit3 Xit 3 = Vit 3 X,t 3 

3 3 



x'i = (3vit* - xd 3)t's + (2a;it"3~ - 3«it')t'3 



/ 2, I 2, / i 1 , 1 2, , 2 



The equations ( A13) and (A14) are the equations of the trajectories. Both coordinates and velocities 
of the particle at a given time t' are expressed through coordinates and velocities at time t ji]] , J24| . 
These equations are utilized in equation ([)]) making it possible to take into account the cosmological 
expansion of the background. 



(A7) 
(A8) 

(A9) 
(A10) 



(All) 
(A12) 

(A13) 

(A14) 
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FIG. 2. The traditional free-streaming considerations. The Jeans length (in dashed-dotted lines); the free-streaming length 
(in solid lines) and the horizon length (in dashed lines) are given as a function of the red shift (qualitative drawing) . Perturba- 
tions corresponding to a wavelength Ai and A3, by crossing the Aj eans line, give origin to a process of collapse, leading to the 
onset of structures. The perturbations with wavelength A2 and A4, by crossing the free-streaming line, are dissipated. On the 
lower right side we consider perturbations originating in the nonrelativistic regime (see text). Xjeansmax and z n . r . depend on 
the fermion mass. 
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FIG. 4. Properties of the instability shell in the nonrelativistic regime. The horizontal lines corresponding to wavelengths Ai 
and A3 and attaining the Jeans length curve (dashed-dotted lines) are the surviving perturbations able to generate structures, 
while A2 and A4 are dissipated by the free-streaming (rising solid lines). At redshift 2 < z n . r . (the vertical line) the instability 
shell (represented by the region with dashed lines) is bounded above by the \j ea ns and below by the solid line defined by 
Aj n / = X ^g s ■ The effects of the free streaming are hampered in the instability shell. The horizon length is indicated by the 
dashed lines. On the lower right site we consider perturbations originating in the nonrelativistic regime (see text). 
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